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Abstract 

We theoretically study the induced odd-frequency pairing states in ballistic normal metal/superconductor 
(N/S) junctions where a superconductor has even-frequency symmetry in the bulk and a normal metal layer 
has an arbitrary length. Using the quasiclassical Green's function formalism, we demonstrate that, quite 
generally, the pair amplitude in the junction has an admixture of an odd-frequency component due to the 
breakdown of translational invariance near the N/S interface where the pair potential acquires spatial de- 
pendence. If a superconductor has even-parity pair potential (spin-singlet s-wave state), the odd-frequency 
pairing component with odd-parity is induced near the N/S interface, while in the case of odd-parity pair 
potential (spin-triplet px-'wawe or spin-singlet dxy-^ave) the odd-frequency component with even-parity is 
generated. We show that in conventional s-wave junctions, the amplitude of the odd-frequency pairing 
state is strongest in the case of fully transparent N/S interface and is enhanced at energies corresponding 
to the peaks in the local density of states (LDOS). In p^- and dxy-^ave junctions, the amplitude of the 
odd-frequency component on the S side of the N/S interface is enhanced at zero energy where the midgap 
Andreev resonant state (MARS) appears due to the sign change of the pair potential. The odd-frequency 
component extends into the N region and exceeds the even-frequency component at energies corresponding 
to the LDOS peak positions, including the MARS. At the edge of the N region the odd-frequency com- 
ponent is non-zero while the even-frequency one vanishes. We show that the concept of odd-frequency 
pairing is a useful tool to interpret a number of phenomena in non-uniform superconducting systems, like 
McMillan-Rowell and midgap Andreev resonance states. 
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I. INTRODUCTION 



Odd-frequency superconducting pairing state, characterized by pair amplitude which is an odd 
function of energy or Matsubara frequency, was first predicted by Berezinskii [1] and has been at- 
tracted a lot of interest recently. Although the existence of odd-frequency pairing in bulk uniform 
systems is not fully established yet Mbi 4, 5, 6], there is a number of proposals to realize it in su- 
perconducting junctions. The realization of the odd-frequency pairing state without finite pair po- 
tential was proposed by Bergeret, Volkov and Efetov in Ref.|7|in ferromagnet/superconductor het- 
erostructures with inhomogeneous magnetization and several related works have been presented 



up to now 



. In particular, it was predicted that the local density of states (LDOS) in the 
ferromagnet is enhanced in the presence of the odd-frequency pairing nlW . 

Recently, it was shown that the odd-frequency pairing state is possible even without magnetic 
ordering. Two of the present authors predicted that the odd-frequency pair amplitude can be in- 
duced in a diffusive normal metal attached to a spin-triplet superconductor [|l2n . According to this 
study, the orig in of the anomalous proximity effect specific to spin-triplet p-wave superconductor 
junctions [|l3ll is the realization of the odd-frequency pairing state in the diffusive normal metal. 



1511 into 



It is also clarified that the penetration of the midgap Andreev resonant state (MARS) [|14l 
the diffusive normal metal is the manifestation of the existence of the odd-frequency spin-triplet 
s-wave superconducting state. The MARS is the well-known resonant state specific to uncon- 
ventional superconductors with sign change of the pair potential on the Fermi surface and was 
observed experimentally in various materials lll6l]. 
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1811 that due to spatial variation of the pair po- 



Furthermore, it was predicted very recently 

o 

tential near a N/S interface [11911 . the odd-frequency pairing state can be induced even in a conven- 
tional ballistic N/S system without spin-triplet ordering. By studying infinite normal metal/infinite 
superconductor (N/S) junctions, it was shown that, quite generally, the spin-singlet even-parity 
(spin-triplet odd-parity) pair potential in the superconductor induces the odd-frequency pairing 
component with spin-singlet odd-parity (spin-triplet even-parity) near the N/S interface [|l7|l . The 
magnitude of the induced odd-frequency component is enhanced in the presence of the MARS due 
to the sign change of the anisotropic pair potential at the interface. In Ref. [3. only the genera- 
tion of the odd-frequency component at the S side of the N/S interface was studied by two of the 
authors. Therefore the question remains how this component extends into the N region and how 
it manifests itself in the properties of the normal metal. In a semi-infinite ballistic normal metal 
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attached to a superconductor, the LDOS normalized by its value in the normal state is always 
unity. This well-known property of LDOS is due to the absence of interference between electrons 
and Andreev reflected holes in a semi-infinite N metal. In this case, the LDOS cannot be used 
to characterize the superconducting correlations in a normal metal. Thus, in order to understand 
manifestations of the induced odd-frequency pairing state in the N metal in a much more clear 
way, it is necessary to study junctions with finite length of the N region. 

In the present paper, using the quasiclassical Green's function theory, we study the pair am- 
plitude and the LDOS at the N/S interface when the N region has finite thickness L. The spatial 
dependence of the pair potential is determined self-consistently. For the convenience of the ac- 
tual numerical calculation, we have used the boundary condition in the Ricatti parametrization 
I 20I1 . The superconductor is assumed to have the conventional even-frequency pairing state in the 
bulk, being in the spin-singlet even-parity state (s-wave or dxy-vjave symmetry) or in the spin- 
triplet odd-parity state (p^-v^ave symmetry). We show that, quite generally, the spatial variation 
of the pair potential and the proximity effect lead to the generation of the odd-frequency com- 
ponent near the N/S interface and on the N side. Moreover, when the superconductor is in the 
even-parity (odd-parity) state, the resulting odd-frequency component is odd-parity (even-parity) 
in order to conserve the spin component. In the absence of the MARS, like in the case of spin- 
singlet s-wave junctions, the magnitude of the odd-frequency component of the pair amplitude is 
suppressed when the transmission coefficient through the interface decreases. The resulting odd- 
frequency pair amplitude has its maximum value at the interface. At the edge of the N region, the 
odd-frequency component is always absent as well as in the S region far away from the interface. 
The LDOS is suppressed around e = 0, where e is the quasiparticle energy measured from the 
Fermi level. For large magnitude of L, the resulting LDOS has an oscillatory e dependence. The 
amplitude of the odd-frequency pair amplitude can exceed that of the even-frequency one at some 
e values. For the case of spin-triplet ^^.-wave and spin-singlet d^y-wave junctions, the amplitude 
of the odd-frequency component at the S side of the N/S interface is much larger than that of the 
even-frequency pair amplitude. This is due to the fact that the presence of the MARS at the inter- 
face enhances the amplitude of the odd-frequency paring state as shown in Ref. |l7. At the 
edge of the N region, the even-frequency component is always absent and only the odd-frequency 
component is nonzero. At e = the resulting odd-frequency component always exceeds the 
even-frequency one. 

The organization of the present paper is as follows. In section 2, we introduce the model and the 
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quasiclassical Green's function formalism. In section 3, the results of the numerical calculations 
are discussed for the case of spin-singlet s-wave, spin-triplet p^^-wave and spin-singlet rf^^j^-wave 
junctions. In section 4, the conclusions and outlook are presented. 



II. MODEL AND FORMULATION 



In the following, we consider a N/S junction as the simplest example of non-uniform supercon- 
ducting system without impurity scattering. Both cases of spin-triplet odd-parity and spin-singlet 
even-parity symmetries are considered in the superconductor. In the spin-triplet superconductor 
we choose Sz = for simplicity. We assume a thin insulating barrier located at the N/S in- 
terface (x = 0) with N (— L < X < 0) and S (x > 0) modeled by a delta function HS{x), 
where H is the strength of the delta function potential. The length of the normal region is L. 
The reflection coefficient of the junction for the quasiparticle for the injection angle 9 is given by 
R = Z^/{Z^ + Acos^e) with Z = 2H/vf, where 6 {-tt /2 < 6 < 7r/2) is measured from the 
normal to the interface and fi;' is the Fermi velocity. 

The quasiclassical Green's functions [21] in a normal metal (N) and a superconductor (S) in 
the Matsubara frequency representation are parameterized as 

^« = fith + f^lh + ^«f3, {g^Sf = i (1) 

where the subscript i{= N, S) refer to N and S, respectively. Here, fj (j = 1, 2, 3) are Pauli 
matrices and i is a unit matrix. The subscript +(— ) denotes the left (right) going quasiparticles 



[|22ll . Functions satisfy the Eilenberger equation [|23ll 

tvpJl^ = T[H±,gt^] (2) 

with 

H± = iUnTs + iA±(x)r2. 

Here vp^ is the x component of the Fermi velocity, Un = 27rT(?2+l/2) is the Matsubara frequency, 
n is an integer number and T is temperature. A+(x) (A_(x)) is the effective pair potential for left 
(right) going quasiparticles. In the N region, A±(x) is set to zero due to the absence of a pairing 
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interaction in the N metal. The above Green's functions can be expressed as 



/« = ±^(F« + - Z}«F«), (3) 

/« = -(F«-D«)/(l-^M^ 
^« = (1 + - D'i^F^^), 

where D^J_\x) and F^\x) satisfy the Ricatti equations [20] in the N region 



vpAD^r (^) = -2a;„Df ^ (x) (4) 
VF^d^Fi''\x) = 2uj^Fi''\x), (5) 

and in the S region, 

VF.d^D^^\x) = -A±{x)[l - {D^f\x))''] + 2uo^Df\x) (6) 
VF^F^ix) = -A±(x)[l - {Fi'\x)f] - 2u;^Fi'\x), (7) 

respectively. 

The boundary conditions at the edge of N region, x = —L, have the form 

i^f = -Di^V^), f':''\-L) = -D'fX-L) (8) 



The boundary conditions at the N/S interface, a; = 0, are 

{l-R)Df\<:)) + [R + L 

[1 + RDf\{])D^!'\<d)] + (1 - R)D'^'\<d)D'^\^) 



p(5).,. ^ (1 - r)d^^\q) + [R + Z^f H0)Di^)(0)]Df (0) 



and 

pjN),^. _ (1 - (0) + \R + (o)DL^Ho)]Df Ho) (10) 

where R is the reflection coefficient at the interface. Since there is no pair potential in the N region, 
the solutions for the spatial dependence of above functions can be easily found 

d'^^\x) = -At-\ Df\x) = -Bt-' 
Ff\x) = At, F[''\x) = Bt, 

with t = exp[(x + L)/^] and ^ = hvFxf^l^nl- The constants A and B are given by 

-2(l-i?)Df (0)to ^^^^ 
A + v/a2 + 4(1 - R)HlD^_^\Q)D^f\0) 
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A + ■^A2 + 4(1 - R)HID^^\())D^^\q) 

with 

A = tl{l + RDf{<d)D'^^\<S)) -{R + D^_^\0)d':^\0)) 
with to = exp(L/^). After simple manipulation, we obtain /^^^\ /2^'' and g^'^ 

A^f = ~z{At - B/t)/{l + AB), = z{Bt - A/t)/{l + AB), 

(13) 

/if = -{At + B/t)/{l + AB), fZ^ = -{Bt + A/t)/{l + AB), 



Note that as follows from Eq. (O, functions gZ"* and (yfi^'' are spatially-independent. 

Here, we consider the situation without mixing of different symmetry channels for the pair 
potential. Then the pair potential A±(x) is expressed by 

A±(x) = A(x)<l>±(^)0(x) (15) 

with the form factor $±(6') given by ^±{6) = 1, ±sin2^^, and ±cos^^ for s-wave, d^^-wave, 
and p3;-wave superconductors, respectively. The pair potential A(x) is determined by the self- 
consistent equation 

A(x) = ^ V / d9G{9)f2+ (16) 



n>l 



n>0 



with G(6') = 1 for s-wave case and G{6) = 2$(6') for other cases, respectively n24\ . Tc is the 
transition temperature of the superconductor. The condition in the bulk is A(oo) = Aq. Since the 
pair potential A(x) is a real quantity, the resulting fi± is an imaginary quantity and f2± is a real 
one. 

Before performing actual numerical calculations, we now discuss general properties of the pair 
amplitude. In the following, we explicitly write f[^]. = flj_{ujn,0), /^^ = fili^niS), F±:'' = 
F^^ {un, 9) and d'^^ = D^^ (c<j„, 9). For the limit x = oo, we obtain 



/iVl^n, 0) = 0, /2± (Wn, 0) = ===. (17) 

V2 + A2$2^(0±) 



Note that flj.{ujn, 9) becomes finite due to the spatial variation of the pair potential and it does 
not exist in the bulk. From Eqs. (U) and Q, we can show that D^J_\—uJn, 9) = 1/D^^\(jjn, 9) and 



F^\—uJn, 0) = l/F^:\uJn, 0). After simple manipulation, we obtain 

0) = -ftk-^n, e), /2k, e) = /«(-c^n, e), 



(18) 



for any x. It is remarkable that functions f[^]_{ujn, 9) and /, 
and even-frequency components of the pair amplitude, respectively |8, 17]. Function 6) 
describes the odd-frequency component of the pair amplitude penetrating from the superconductor. 

Next, we discuss the parity of these pair amplitudes. The even-parity (odd-parity) pair 
amplitude should satisfy the following relation /j2(to'„,6') = /j!)(u;„, — 6^) 6*) = 



correspond to odd-frequency 



-6)], with j = 1, 2. For an even-parity (odd-parity) superconductor, $±(—6*) = ^^{0) 
— $zp(6')]. Then, we can show that for the even-parity case 



{i), 



F^\9) 



(19) 



and for the odd-parity case 



D 



respectively. 

The resulting 6) and /2±(cu„ 



D^\9), F 



F^\9) 



satisfy 

/(2(c.„,^) = -/« 



(20) 



for an even-parity superconductor and 



/: 
J'. 



(21) 



for an odd-parity superconductor, respectively [|l7h . Note that the parity of the odd-frequency 



component /2(c<j„, 9) is always different from that in the bulk superconductor. 

As shown above, the odd-frequency component f[^^{ujn,9) is purely an imaginary quantity. 
The underlying physics behind this formal property is the follows. Due to the breakdown of trans- 
lational invariance near the N/S interface, the pair potential A(x) acquires a spatial dependence 
which leads to the coupling between even-parity and odd-parity states. Since the bulk pair poten- 
tial has an even-frequency symmetry, the Fermi-Dirac statistics requires that the order parameter 
component induced near the interface should be odd in frequency. The phase of the induced pair 
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amplitude undergoes a 7r/2 shift from that in the bulk S thus removing internal phase shift between 
the even- and odd-frequency components and making the interface-induced state compatible with 
the time reversal invariance. As a result, function f^^i^n, 0) becomes a purely imaginary quantity 



Let us now focus on the values of the pair amplitudes at the edge of N region (at x = —L). We 
concentrate on two extreme cases with (I) ^+{6) = ^-{9) and (II) ^+{9) = —^-{9). In the case 
(I), the MARS is absent since there is no sign change of the pair potential felt by the quasiparticle 
at the interface. Then the relation D^_^^ = holds. On the other hand, in the case (II), the 
MARS is generated near the interface due to the sign change of the pair potential and the relation 



D^_^^ = —D^_^^ is satisfied [|14l1 . At the edge x = —L, it is easy to show that F^'' = —D]^''' for 



(JV) _ n(7V) 



the former case and F^'^ = for the latter one. As a result, f[^^ = for the case (I) and 
/^(f ) = for the case (II), respectively. Thus we can conclude that in the absence of the MARS 
only the even-frequency pairing component exists at x = —L, while in the presence of the MARS 
only the odd-frequency one. 

In order to understand the angular dependence of the pair amplitude in a more detail, we de- 
fine /P and ff^ for -7r/2 < 9 < 37r/2 with /jg^ = /|(^)+(^) for -7r/2 < 9 < 7i/2 and 
A^(2) ~ /i(2)-(^ — 6*) for 7r/2 < < 37r/2. We decompose /|^2) ii^to various angular momentum 
component as follows, 

fm = E S'^^"^' + E (22) 

m m 

with m = 2/ + 1 for odd-parity case and m = 21 for even-parity case with integer / > 0, where / is 
the quantum number of the angular momentum. Here, CnS"^^^ and Sm^'^^^ are defined for all x. It is 
straightforward to show that the only nonzero components are (1) and c!^]]^^ for even-parity 

(2) 

superconductor without sign change at the interface (i.e., s-wave or (ij.2_y2-wave), (2) S2i_^_2 and 
5*2^1]^ for (ij.y-wave, (3) Cg^l^ and 02]^ for j9^.- wave, and (4) 5*2^11 and 82]^ for p^^-wave junctions, 
respectively. The allowed angular momenta for odd-frequency components are 2/ + 1, 2/ + 1, 2/, 
and 2/ + 2 corresponding to each of the above four cases. 

In order to get better insight into the spectral property of the odd-frequency pair amplitude, we 
perform an analytical continuation from the Matsubara frequency cOn to the quasiparticle energy 
e measured from the chemical potential. The retarded Green's function corresponding to Eq. (1) 
is defined as = fifn + /g^f^ + g^^%. One can show that f^'^i-e) = -[/Jf (5)]*, 
f^l^{-6) = [/2±^ (£)]*, and ^^^^(-5) = [g^±^^{6)]*. TheLDOSp(£) at the N/S interface at x = 



normalized to its value in the normal state is given by 



p(e) = / d^Real( 



J-w/2 



/.7r/2 




(23) 



In the following, we self-consistently calculate the spatial dependence of the pair potential and 
the pair amplitude in the Matsubara representation. After that we calculate the spectral properties 
of pair amplitudes and LDOS. For actual calculations, we choose spin-singlet s-wave, spin-triplet 
p^.-wave and spin-singlet rf^^-wave state in a superconductor and fix temperature T = O.OSTc. 
The length of the normal region L is measured in units of Lq = vp/2ttTc. 

III. RESULTS 

A. s-wave pair potential 

First we focus on the s-wave superconductor junctions as shown in Fig. [T] By changing the 
length L of the N region and the transparency at the interface, we calculate the spatial dependence 
of the pair potential and the pair amplitudes in the Matsubara frequency representation. We only 
concentrate on the lowest angular momentum of the even-frequency pair amplitude Cq . As re- 
gards the odd-frequency pair amplitudes, we focus on the C[^\ C^^^ and C^'^^ components which 
all have odd-parity and depend on 9 as cos 9, cos 36' and cos 59, respectively, and correspond to 
p^-wave, /i-wave and /ii-wave components shown in Fig. 1. In all cases, even-frequency compo- 
nent is constant in the S region far away from the interface and the corresponding odd-frequency 
components are absent. The s-wave pair potential is suppressed for the fully transparent case 
(Z = 0), while it is almost constant for low transparent case (Z = 5). It does not penetrate into 
the N region due to the absence of the attractive interaction in the N metal. On the other hand, 
in all considered cases the spatial variation of the even-frequency s-wave pair amplitude is rather 
weak in the S region, while in the N region it is strong for Z = and is reduced for Z = 5 since 
the proximity effect is weaker in the latter case. The odd-frequency component always vanishes 
at a; = — L and does not have a jump at the N/S interface even for nonzero Z. Its amplitude is 
strongly enhanced near the N/S interface especially for fully transparent junctions. Note that not 
only the ^^.-wave but also /i-wave and /ii-wave have sufficiently large magnitudes as shown in 
Figs.dja) and [2c). With the decrease of the transparency of the N/S interface, the odd-frequency 
components are suppressed as shown in Figs.IIJb) and[2d). 
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FIG. 1 : (Color online) Spatial dependence of the normalized pair potential, even-frequency pair amplitude 
and odd-frequency components of the pair amplitude for s-wave superconductor junctions. Here, we choose 
^ = tip'/Ao in the S region (x > 0) and ^ = Lq = vf/2itTc in the N region. The pair amplitudes 
Cq^\ c[^\ C^^\ and C^^^ are denoted as even s-wave, odd px-wave, odd /i-wave, and odd /ii-wave pair 
amphtudes. (a) Z = 0, L = Lq, (b) Z = 5, L = Lq, (c) Z = 0, L = SLq, and (d) Z = 5, L = 5Lq, 
respectively. 

In order to understand the proximity effect in more detail, we look at the resulting LDOS and 
the spectral properties of pair amplitudes in real energy e. We focus on the even-frequency s-wave 
component of the pair amplitude Cq^^ and on the odd-frequency pa;-wave pair amplitude C[^^ on 
the S side of the N/S interface and on the edge of the N region. As follows from Eq. ([14)) . the 
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FIG. 2: (Color online) Energy dependence of the LDOS and the pair amplitudes in s-wave junctions with 
L = Lq and Z = 0. (a) The LDOS normalized by its value in the normal state. The solid line: LDOS 
on the S side of the N/S interface, the dotted line: LDOS in the N region. Energy dependence of the real 
(solid line) and imaginary (dotted line) part of (b) even-frequency s-wave pair amplitude on the S side of 
the N/S interface, (c) odd-frequency pa;-wave pair amplitude on the S side of the N/S interface and (d) 
even-frequency s-wave pair amplitude at the edge of the N region. 

LDOS is independent of the coordinate x in the N. For Z = and L = Lq, the LDOS has a 
V-shaped structure. There is no jump of the value of the LDOS at the N/S interface. The even- 
frequency pair amplitude at the N/S boundary on the S side is shown in Fig. I2tb). Its real part is 
an even function of e while its imaginary part is an odd function of e. The corresponding odd- 
frequency one is plotted in Fig. Htc). In contrast to the even-frequency component [Fig. Htb)], the 
real (imaginary) part of the pair amplitude is an odd (even) function of e. The pair amplitude is 
enhanced around e ~ ±0.6Ao where the LDOS have a peak. At the N/S boundary (x = —L), only 
the even-frequency component exists. The line shape of the pair amplitude shown in Fig. WiA) is 
similar to that in Fig. |2lb). 

In Fig. [3l the corresponding plots for L = Lq and Z = 5 are shown. The LDOS on the S 
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FIG. 3: (Color online) Same as Fig.|2j but with L = Lq and Z = 5. 

side of the N/S interface has a U-shaped DOS similar to bulk DOS. On the other hand, in the N 
region, the LDOS has a different value due to the discontinuity at the N/S interface. The LDOS in 
N also exhibits the minigap structure which scales with the interface transparency, in accordance 
with the well-known McMillan model of proximity effect in conventional superconducting N/S 
junctions [|25l] . The magnitude of the real part of the even-frequency component on the S side of 
the N/S interface exceeds the magnitude of the imaginary part as seen from Fig.[3];b) for < Aq. 
The magnitude of the odd-frequency part is small as compared to that of the even-frequency one. 
As seen from Fig.[3tc), the real part of the odd-frequency component has a minigap structure and 
the imaginary part has a dip and peak structure in contrast to the case of the even-frequency one 
[see Fig.[3tb)]. At the N/S boundary (x = —L), only the even-frequency component exists. The 
real part of the even-frequency component at x = —L has a peak around e = [see Fig. (Sd)]. 
The width of this peak is of the same order as that of the dip of LDOS. As compared to the 
corresponding case of Z = 0, the proximity effect in the N region is only essential at low energy 
e. 

It is also interesting to consider the case of large width of the N region. Here, we concentrate 
on the situation when the N/S/ interface is fully transparent (Z = 0) and L = 5Lq. In this case the 
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FIG. 4: (Color online) Same as Fig.[2l but with L = 5Lq and Z = 0. 



LDOS in the N region and at the N/S interface coincide with each other as seen from Fig.lH The 
LDOS has multiple peaks due to the existence of the multi-sub gap structures due to electron-hole 
interference effects in the N region [|26n . The amplitudes of the corresponding even-frequency and 
odd-frequency components are enhanced at energies e corresponding to the LDOS peak positions, 
while the ratio of this components depends on energy and location in the N region. To clarify this 
point much more clearly, we concentrate on the ratio of the odd- and even-frequency components 
in the N region. According to Eq. (fT3l) . the ratio of the magnitude of the odd-frequency component 



fi^\s, 9) to the even-frequency one f2+'{^, is 



\Ut-t\ 
\l/t + t\ 



2^ 

tan ( (L + 

VFx 



(24) 



At the edge of the N region, x = —L, the odd-frequency component vanishes at all energies. 
On the other hand, very interesting situation occurs at the N/S interface, x = as will be shown 
below. In Fig. [51 we plot this ratio for 9 = and x = 0. It is remarkable that at some energies the 
amplitude of the odd-frequency pair amplitude exceeds that of the even-frequency one. 

Let us clarify the relation between the positions of the bound states and the above ratio of the 
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FIG. 5: Ratio of the pair amplitudes f[^\e,9)/f2^\e,9) on the N-side of the N/S interface in s-wave 
junction as a function of energy e for ^ = and L = BLq. 



odd-to-even pair amplitude. In the limit L >> Lq the bound states are determined by simple 
relation [|26|] 



2L 



(n + 1/2), n = 0,1,2, 



(25) 



Very dramatic situation occurs at the N/S interface, a; = 0: combining the above two equations, we 
obtain that at the LDOS peak positions e = En the ratio of the odd-to-even pair amplitude diverges 



2+ \ 



|tan (7r/2 + 7rn)| 



oo. 



(26) 



That means that at the subgap peak energies the odd-frequency component dominates over 
the even-frequency one at the N/S interface. This is a remarkable property of the odd-frequency 
pairing, which makes it relevant to the classical McMillan-Rowell oscillations in the N/S geometry 
ll26ll . To summarize, we have shown that the odd-frequency component is present even in the 
standard case of a ballistic N/S system, and it dominates at energies when the LDOS has subgap 
peaks. 



B. px-wave pair potential 

Next, we focus on the p^^-wave superconductor junctions as shown in Fig. [6l Similar to the 
case of s-wave junctions, by changing the length of the normal region L and the transparency at 
the interface, we calculate the spatial dependence of the pair potential and the pair amplitudes in 
the Matsubara frequency representation. We only concentrate on the lowest angular momentum 
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of the even-frequency pair amplitude C{ . As regards the odd-frequency pair amplitude, we focus 
on the Cq^\ Cg^'* and C^^^ components where the parity of the odd-frequency components is even. 
These functions correspond to s-wave, (i^2_j^2-wave and (7- wave components in Fig. [6l where 
the 9 dependencies are given by 1, cos26', and cos46', respectively. In all cases, even-frequency 
component is constant in the S region far away from the interface and the corresponding odd- 
frequency components are absent. The p^-wave pair potential is reduced at the N/S interface in 
all cases. For Z = 5, the reduction is significant and the resulting magnitude of the p^-wave 
pair potential is almost zero at the N/S interface. It does not penetrate into the N region due 
to the absence of the attractive interaction in the N metal. The amplitude of the p^^-wave even- 
frequency pair amplitude is reduced towards the N/S interface and monotonically decreases in 
the N region. It does not have a jump at the N/S interface even for nonzero Z and vanishes at 
the edge of the N region (x = —L). On the other hand, the odd-frequency component is always 
nonzero at x = —L and has a jump at the N/S interface for nonzero Z [see Figs.[6tb) and[6td)]. 
The amplitude of the odd-frequency component is strongly enhanced near the S-side of the N/S 
interface. This enhancement is much more significant for the low transparent interface with large Z 
[see Figs.[6tb) and[6td)]. However, for the presently chosen Matsubara frequency ujn = O.OSvrTc 
it cannot penetrate into the N region. On the other hand, for Z = 0, the odd-frequency component 
significantly extends into the N region. Note that not only the s-wave but also d.j.2_y2-wave and 
(yf-wave components have sufficiently large magnitudes as shown in Figs.^a) and[6tc). These pair 
amplitudes are almost constant in the N region. 

In order to get better insight into the spectral property of the odd-frequency pair amplitude, we 

calculate the LDOS and the pair amplitudes as functions of real energy e. We focus on the even- 

(2) 

frequency pj;-wave component of the pair amplitude CI and odd-frequency s-wave component 
of the pair amplitude Cq^^ at the N/S interface on the S side and the N boundary. In the N region, 
the LDOS is independent of x as shown by Eq. ([141) . For Z = and L = Lq [see Fig.|7ta)], the 
LDOS has a zero energy peak (ZEP) due to the formation of the MARS. There is no jump of the 
LDOS at the interface since the interface is fully transparent. The even-frequency pair amplitude 
at S-side of the N/S boundary is shown in Fig.|7tb). Both the real and imaginary parts do not vary 
strongly around e ~ 0. 

Similar to the case of s-wave junctions, the real part of the even-frequency component is an 
even function of e while the imaginary part is an odd function of e. The corresponding odd- 
frequency component is plotted in Fig.|7tc). The real (imaginary) part of the pair amplitude is odd 
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FIG. 6: (Color online) Spatial dependence of the normalized pair potential, even-frequency and odd- 
frequency pair amplitudes for p^-wave superconductor junctions. Here, we choose = vp/ Aq in the S 
region {x > 0) and ^ = Lq = vp/^irTc in the N region. The pair amplitudes c[^\ Cq^\ Cg^'', and C^^ 
are denoted as even p^-wave, odd s-wave, odd d^2_y2-v/ave, and odd g-wave pair amplitudes, (a) Z = 0, 
L = Lq, (b) Z = 5, L = Lq, (c) Z = 0, L = 5Lq, and (d) Z = 5, L = 5Lq, respectively. 

(even) function of e. The amplitude of the odd-frequency pair amplitude is enhanced around e = 
where the LDOS has the ZEP. At the edge of the N (x = —L), only the odd-frequency component 
exists. The line shape of the pair amplitude shown in Fig.Hd) is similar to that shown in Fig.Hc). 

In Fig. [HI the results of corresponding calculation with L = Lq and Z = 5 are shown. Both the 
LDOS at the N/S interface and the edge of the N have a ZEP. In the N region, the LDOS is almost 
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FIG. 7: (Color online) Energy dependence of the LDOS and the pair amplitudes in ^^.-wave junctions with 
L = Lq and Z = 0. (a) The LDOS normalized by its value in the normal state. The solid line: LDOS 
on the S side of the N/S interface, the dotted line: LDOS in the N region. Energy dependence of the real 
(solid line) and the imaginary (dotted line) part of (b) even-frequency px-wave pair amplitude on the S side 
of the N/S interface, (c) odd-frequency s-wave pair amplitude on the S side of the N/S interface and (d) 
odd-frequency s-wave pair amplitude at the edge of the N region. 

unity due to the absence of the proximity effect for \e\ > 0.24Ao [see dotted line in Fig.[8ta)]. The 
LDOS has a ZEP and small peak at e = 0.24Ao. The corresponding real and imaginary parts of 
the even-frequency pair amplitude at the N/S interface also have peaks at this energy (Fig.[8tb)]). 
The amplitude of the odd-frequency component is enhanced as compared to the corresponding 
even-frequency one as shown in Fig. [Htc). At the edge of the N region (x = —L), the amplitude 
of the odd-frequency component is almost zero for \ £ \> 0.24Ao. However, around zero energy, 
the amplitude of the odd-frequency component is drastically enhanced as in the case of S-side of 
the N/S boundary. The penetration of the odd-frequency component occurs only at low energies. 

For the longer normal region with L = 5Lo, the resulting LDOS has the ZEP and a number 
of peaks at finite energies e [see Fig.lHa)]. The even-frequency component of the pair amplitude 
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FIG. 8: (Color online) Same as Fig.|7l but with L = Lq and Z = 5. 

on the S side of the N/S boundary also has multiple peaks. The corresponding odd-frequency 
component has many peaks with amplitudes strongly enhanced around e = 0. Around zero energy, 
the amplitude of the odd-frequency component is much larger than that of the even-frequency 
one [see Fig. IH^c)]. At the edge of the N region, the resulting odd-frequency component has a 
significant amplitude as shown in Fig. Hfd). Finally, we focus on the ratio of the odd- and even- 
frequency components of the pair amplitude, f[^\e, 0)/f2+\^, G)- In Fig. [lOl we plot this ratio 
for = and x = 0. Remarkably, at some energies the odd-frequency pair amplitude exceeds 
that of the even-frequency one. In contrast to the s-wave case, there is a huge peak at e = 
corresponding to the existence of the MARS. 

To summarize, we have shown that when the LDOS has a ZEP, the resulting odd-frequency 
component is enhanced around e, its imaginary part having a ZER It is evident that the odd- 
frequency pairing state is indispensable to understand the proximity effect in ^^.-wave supercon- 
ductor system. 
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FIG. 10: Ratio of the pair amplitudes /}_|_ {e,9)/ (e, 6) as a function of e for L = 5Lq at the N-side of 
the N/S interface for px-v^ave junctions for ^ = 0. 

C. dxy-vfa\e pair potential 

Finally we focus on the d^y-wsve junctions as shown in Fig. \TT\ Similar to the above two 
cases, by changing the length of the normal region L and the transparency at the interface, we 
calculate the spatial dependence of the pair potential and the pair amplitudes in the Matsubara 
frequency representation. Here we only concentrate on the lowest angular momentum of the even- 
frequency pair amplitude s!^\ The s-wave component of the pairing amplitude is absent due 
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to the sign change of the pair potential with respect to the exchange 6 by —6. As regards the 
odd-frequency pair amplitude, we focus on the S^^ and S^^ components where the spatial 
parity of the odd-frequency components is odd. These cases correspond to py-wave, /2-wave 
and /i2-wave components in Fig. [TT] where the 9 dependence is given by sin 9, sin 36', and sin 56', 
respectively. In all cases, the even-frequency component is constant in the S region far away from 
the interface and the corresponding odd-frequency components are absent. The dxy-wave pair 
potential is suppressed at the N/S interface in all cases. For Z = 5, the reduction is significant 
and it is almost zero at the N/S interface. The even-frequency dxy-wa\e pair amplitude is reduced 
towards the N/S interface and monotonically decreases in the N region similar to the case of 
Pa;-wave one [see Figs. [6ta) and [TTT a)1 . It does not have a jump at the N/S interface even for 
nonzero Z and vanishes at the edge of N region (x = —L). On the other hand, the odd-frequency 
component is always nonzero at x = —L and has a jump at the N/S interface. The amplitude 
of the odd-frequency component is strongly enhanced near the S-side of the N/S interface. This 
enhancement is much more significant for the low transparent interface with large magnitude of 
Z [see Figs. fnTb) andHUd)]. On the other hand, for Z = 0, the odd-frequency components 
significantly penetrate into the N region. Note that not only the p^^-wave but also /2-wave and /i2- 
wave components have the sufficiently large magnitudes as shown in Figs. [TIT a) and [TTVc). The 
above pair amplitudes are almost constant in the N region. 

In order to get better insight into the spectral property of the odd-frequency pair amplitude, 
we calculate the LDOS and the pair amplitudes as functions of the real energy e. We focus on 
the even-frequency c/^j^-wave component of the pair amplitude 5*2 and odd-frequency p.y-wave 
component of the pair amplitude S']^^ at the S side of the N/S interface and at the edge of the N 
region, x = —L. The resulting LDOS has the ZEP due to the formation of the MARS. Similar to 
the previously considered cases of s-wave and p^^-wave junctions, the LDOS is independent of x as 
follows from eq. (fT4l) . Here, we choose Z = and L = Lq [see Fig.[T2ta)]. Similar to the p^-'wawe. 
case, the even-frequency pair amplitude at the N/S boundary on the S-side is shown in Fig.[l2tb). 
Both the real and the imaginary parts do not vary strongly around £ ~ 0. The real part of the even- 
frequency component is an even function of e while its imaginary part is an odd function of e. The 
corresponding odd-frequency component is plotted in Fig. [T2l' c). In contrast to the even-frequency 
component [Fig. [I2tb)], the real (imaginary) part of the pair amplitude is an odd (even) function 
of e. The amplitude of the pair amplitude is enhanced around £ ~ where the LDOS has the ZEP. 
At the N/S boundary (x = —L), only the odd-frequency component exists. The line shape of the 
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FIG. 11: (Color online) Spatial dependence of the normalized pair potential, even-frequency and odd- 
frequency pair amplitudes for ^^.j^-wave superconductor junctions. Here, we choose ^ = fir/Ao in the S 
region (x > 0) and ^ = Lq = vf/'^t^Tc in the N region. The pair amplitudes s'^\ and S^^ 

are denoted as even d^^-wave, odd pj^-wave, odd /2-wave, and odd /12-wave pair amplitudes, (a) Z = 0, 
L = Lq, (b) Z = 5, L = Lq, (c) Z = 0, L = 5Lq, and (d) Z = 5, L = 5Lq, respectively. 

pair amplitude as shown in Fig. [721 (1) is qualitatively similar to that in Fig.[T2tc). This qualitative 
behavior of the line shapes is very similar to that for the corresponding p^-v^awe case. 

Similar to the Pa;-wave junction case, when the LDOS has the ZEP, the resulting odd-frequency 
component is enhanced around e. It is evident that the odd-frequency pairing state is indispensable 
to understand the proximity effect in (i^^-wave superconductor system. 
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FIG. 12: (Color online) Energy dependence of the LDOS and the pair amplitudes in dxy-v/ave junctions 
with L = Lq and Z = 0. (a) The LDOS normalized by its value in the normal state. The solid line: LDOS 
on the S side of the N/S interface, the dotted line: LDOS in the N region. Energy dependence of the real 
(solid line) and the imaginary (dotted line) part of (b) even-frequency dxy-'wa.ve pair amplitude on the S side 
of the N/S interface, (c) odd-frequency pj^-wave pair amplitude on the S side of the N/S interface and (d) 
odd-frequency p^-wave pair amplitude at the edge of the N region 

In the following, we comment on the differences between odd-frequency pair amplitudes in 
d^y-wave and p^-wawe junctions. In both cases, the magnitude of the odd-frequency component is 
enhanced at the interface and in the normal region. However, the odd-frequency odd-parity state 
is generated for dxy-vjawe case, while the odd-frequency even-parity state is generated for p^-wave 
case. The s-wave isotropic component which is robust against the impurity scattering [I12I] appears 
only in the latter case. Then the dxy-vjave pair amplitude cannot penetrate into diffusive normal 
metal while the p^-wave one can. Thus we can naturally understand the presence of proximity 



effect with the MARS in p^,- wave junctions [I12|,[l31] and its absence in rfa-^-wave junctions 1112 . 
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It is instructive to relate the LDOS anomalies in dxy-wave and pa;-wave junctions to the magni- 
tude of the odd-frequency pairing component. According to Eq. (fT3l) . in d^y- and p^j-wave junc- 
tions the ratio of the magnitude of the odd-frequency component fi^\e, 9) to the even-frequency 
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onef^^\e,0) is 



\l/t + t\ 
\l/t-t\ 



cotan (L + x 

\VFx 



(27) 



It follows from the above expression that at the edge of the N region, x = —L, the odd- 
frequency component dominates at all energies. On the other hand, at the N/S interface, x = 0, 
the odd-frequency component dominates at energies e = En corresponding to the LDOS has peak 
positions 

T^VFxn 

£n= ^ n = 0,1,2,... (28) 



2L 



For n = Eq. (1281) describes the mid-gap Andreev bound state and higher n correspond to the 
subgap resonances for large N region thickness. Therefore, we can conclude that in (ia-^-wave and 
p^-wave junctions the odd-frequency component dominates over the even-frequency one at the 
N/S interface at the energies corresponding to the LDOS peak positions, including the prominent 
zero-energy peak (MARS). Moreover, the odd-frequency component always dominates at the edge 
of the N region, x = —L, where the breaking of translational invariance is the strongest because 
of sign change of the pair amplitude at that point. 



IV. CONCLUSIONS 



In summary, using the quasiclassical Green's function formalism, we have shown that the odd- 
frequency pairing state is ubiquitously generated in the normal metal/superconductor (N/S) ballis- 
tic junction system, where the length of the normal region is finite. It is shown that the even-parity 
(odd-parity) pair potential in the superconductor induces the odd-frequency pairing component 
with spin-singlet odd-parity (spin-triplet even-parity). As regards the symmetry of the supercon- 
ductor, we have chosen typical three cases, spin-singlet s-wave, spin-triplet p^^-wave and spin- 
singlet dxy-wave. In the latter two cases, mid gap Andreev resonant state (MARS) appears at 
the N/S interface. Even for conventional s-wave junctions, the amplitude of the odd-frequency 
pairing state is enhanced at the N/S interface with fully transparent barrier. By analyzing the spec- 
tral properties of the pair amplitudes, we found that the magnitude of the resulting odd-frequency 
component at the interface can exceed that of the even-frequency one. For the case of p^-wave 
and (ia;y-wave junctions, the magnitude of the odd-frequency component at the S side of the N/S 
interface is significantly enhanced. The magnitude of the induced odd-frequency component is 
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enhanced in the presence of the midgap Andreev resonant state due to the sign change of the 
anisotropic pair potential at the interface. The LDOS has a zero energy peak (ZEP) both at the 
interface and in the N region. At the edge of the N region, only the odd-frequency component is 
non-zero. 

The underlying physics behind these phenomena is related to the breakdown of translational 
invariance near the N/S interface where the pair potential A(x) acquires a spatial dependence. As a 
result, an odd-frequency component is quite generally induced near the interface. The breakdown 
of translational invariance is the strongest when the pair potential changes sign upon reflection like 
in the case of p^-wave and d^j^-wave junctions, then the magnitude of odd-frequency component is 
the largest. Moreover, the phase of the interface-induced odd-frequency component has a 7r/2 shift 
from that in the bulk of S. Therefore, as shown above, the odd-frequency component /j^2(^n) ^) 
becomes purely imaginary quantity and the peak structure in the LDOS naturally follows from the 
normalization condition. 

We have also shown that in the N/S junctions with s-wave superconductors the classical 
McMillan-Rowell oscillations 12611 can also be reinterpreted in terms of odd-frequency pairing. 
As follows from Eq.(l26l). at the energies corresponding to the subgap peaks in the N/S junction, 
the odd-frequency component dominates over the even-frequency one. This is remarkable appli- 
cation of the odd-frequency pairing concept when one can re-interpret the well-known resonance 
structure. 

In the present study, we clarified the details of proximity effect of the odd-frequency pairing 
state induced at the N/S boundary. We reinterpreted the appearance of the MARS in terms of the 
enhanced odd-frequency pair amplitude. Though we explicitly studied the N/S junctions only, the 
odd-frequency pairing state is also expected near impurities and within Abrikosov vortex cores, 
where the amplitude of the pair potential is reduced. The present result indicates the ubiquitous 
presence of odd-frequency pairing states because most of real superconductors are not uniform. 
That means that the odd-frequency pairing is not at all a rare situation as was previously assumed. 
Thus we believe that the odd-frequency pairing may become an important concept in understand- 
ing the physics of non-uniform superconducting systems. 

In the present paper, the proximity effect is studied in the ballistic limit. In the present case, the 
enhanced odd-frequency pair amplitude appears in the N region both for pa;-wave and (i^j^-wave 
junctions. It is very interesting to study in the intermediate regime JlS] since the parity of these 
states are different. In the diffusive limit, the proximity effect survives only for the former case. 
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There are several interesting phase coherent effects relevant to MARS Hioi 3oll . These preexisting 
phenomena can be reinterpreted in terms of the odd-frequency paring state. 
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